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Abstract 

A direct reconstruction algorithm for complex conductivities in W'^'°°{Q), 
where is a bounded, simply connected Lipschitz domain in R'^, is pre- 
sented. The framework is based on the uniqueness proof by Francini [In- 
verse Problems 20 2000], but equations relating the Dirichlet-to-Neumann 
to the scattering transform and the exponentially growing solutions are 
not present in that work, and are derived here. The algorithm constitutes 
the first D-bar method for the reconstruction of conductivities and per- 
mittivities in two dimensions. Reconstructions of numerically simulated 
chest phantoms with discontinuities at the organ boundaries are included. 

1 Introduction 

The reconstruction of admittivies 7 from electrical boundary measurements is 
known as the inverse admittivity problem. The unknown admittivity appears as 
a complex coefficient 7(2;) = (t(z) + iuje{z) in the generalized Laplace equation 

V-(7(z)Vu(z)) = 0, zen, = (1) 

where u is the electric potential, a is the conductivity of the medium, e is the 
permittivity, and uj is the temporal angular frequency of the applied electromag- 
netic wave. The data is the Dirichlet-to-Neumann, or voltage-to-current density 
map defined by 

du 
ov 



A^u = 7- 



(2) 

an 



where u G H^ip) is the solution to (UJ. By the trace theorem A^ : H^/^{diVj 
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In this work we present a direct reconstruction algorithm for the admittiv- 
ity 7. The majority of the theory is based on the 2000 paper by Francini [20] in 
which it is estabhshed that if cr, e € W^'°°{^1), where $7 is a bounded domain in 

with Lipschitz boundary, then the real-valued functions a and e are uniquely 
determined by the Dirichlet-to-Neumann map, provided that the imaginary part 
of the admittivity is sufficiently small. The proof in is based on the D-bar 
method and is nearly constructive, but equations linking the scattering trans- 
form and the exponentially growing solutions to the Dirichlet-to-Neumann data 
are not used in the proof, and so it does not contain a complete set of equa- 
tions for reconstructing the admittivity. In this work, we derive the necessary 
equations for a direct, nonlinear reconstruction algorithm for the admittivity 7. 
Furthermore, we establish the existence of exponentially growing solutions to 
([1]), which prove to be useful in relating the Dirichlet-to-Neumann data to the 
scattering transform. The reconstruction formula in |20| is for the potential Q^, 
whose relationship to 7 is described below. We provide a direct formula for 7 
from the D-bar equations in 20 , which is computationally advantageous as well. 

The inverse admittivity problem has an important application known as 
electrical impedance tomography (FIT). The fact that the electrical conductivity 
and permittivity vary in the different tissues and organs in the body allows one 
to form an image from the reconstructed admittivity distribution. In the 2-D 
geometry, FIT is clinically useful for chest imaging. Conductivity images have 
been used for monitoring pulmonary perfusion [51 1221 147| , determining regional 
ventilation in the lungs [23l [211 HI] j and the detection of pneumothorax [15] , for 
example. In three dimensions, conductivity images have been used, for instance, 
in head imaging [50l [49] and knowledge of the admittivity has been applied to 
breast cancer detection [H |30l |3l] . 

Reconstruction algorithms based on a least-squares approach that recon- 
struct permittivity include [HI [H [29] . The aforementioned algorithms are iter- 
ative, whereas the work presented here is a direct method that makes use of ex- 
ponentially growing solutions, or complex geometrical optics (CGO) solutions, 
to the admittivity equation. The steps of the algorithm are to compute these 
CGO solutions from knowledge of the Dirichlet-to-Neumann map, to compute a 
scattering transform matrix, to solve two systems of B (D-bar) equations in the 
complex frequency variable k for the CGO solutions to a related elliptic system, 
and finally to reconstruct the admittivity from the values of these solutions at 
fc = 0. In this work, we provide a complete implementation of this algorithm 
and present reconstructions of several numerical phantoms relevant to medical 
FIT imaging. The phantoms we consider here are discontinuous at the organ 
boundaries, which is actually outside the theory of the algorithm. The work [25] 
contains computations of smooth admittivities and validates our formulas and 
computations by comparing the results of the intermediate functions (CGO so- 
lutions and scattering transforms) with those computed from knowledge of the 
admittivity. 

We briefly review the history of results using CGO solutions on the inverse 
conductivity problem in dimension 2. The inverse conductivity problem was first 
introduced by A. P. Calderon [12] in 1980, where he proved that, in a linearized 
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version of the problem, the Dirichlet-to-Neumann map uniquely determines the 
conductivity, and he proposed a direct reconstruction method for this case. 
An implementation in dimension two for experimental data is found in In 
1996, Nachman [33] presented a constructive proof of global uniqueness for 
twice differentiable conductivities using D-bar methods. The D-bar algorithm 
following from [331 131] has been applied to simulated data in [35] 1301 HSl I3Z] and 
to experimental data on tanks and in vivo human data in [27l [28l HT] [17] . While 
the initial scattering transform was regularized using a Born approximation, a 
more recent paper [39] contains a full nonlinear regularization analysis, including 
estimates on speed of convergence in Banach spaces, for twice differentiable 
conductivities. The regularity conditions on the conductivity were relaxed to 
once-differentiable in [10]. The proof uses D-bar techniques and formulates 
the problem as a first-order elliptic system. A reconstruction method based 
on [TU] can be found in [311 1331 133] ■ Francini [3^ provided a proof of unique 
identifiability for the inverse admittivity problem for cr, e G W'^'°°{^1), with uj 
small. Her work provides a nearly constructive proof based on D-bar methods 
on a first-order elliptic system similar to that in [TU] . A non-constructive proof 
that applies to complex admittivities with no smallness assumption is found in 
[TTj . Astala and Paivarinta provide a CGO-based constructive proof for real 
conductivities a S L°°(il), and numerical results related to this work can be 
found in [113] . 

The paper is organized as follows. In Section [2] we describe the direct re- 
construction algorithm, which is comprised of boundary integral equations for 
the exponentially growing solutions to (jlj involving the Dirichlet-to-Neumann 
data, boundary integral equations relating those CGO solutions and the CGO 
solutions ^E" of the first order system, equations for the scattering transform in- 
volving only the traces of \E', the Bk equations established in [20], and the direct 
reconstruction formula for Qj and thus 7. Derivations of the novel equations are 
found in this section. Section [3] describes the numerical implementation of the 
algorithm. Results on noisy and non-noisy simulated data of a cross-sectional 
chest with discontinuous organ boundaries are found in Section |31 

2 The Direct Reconstruction Algorithm 

In this section we will provide the equations for the direct reconstruction algo- 
rithm, completing the steps for the proof in [2D| to be completely constructive. 
In particular, boundary integral equations relating the CGO solutions to the 
Dirichlet-to-Neumann (DN) map are derived. 

Let C be a bounded open domain with a Lipschitz boundary. Through- 
out we assume that there exist positive constants ctq and /3 such that 



a{z) > (To, z e 1^ C M' 



2 



(3) 



and 



lkllwi.~(n), l|e||M/i.~(n) < P- 



(4) 
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We extend a and e from 17 to all of such that ct = 1 and e = outside a 
ball with fixed radius that contains $7, and ([3|) and (|4]) hold for all of . In 
fact, all that is required is that 7 is constant outside that ball of fixed radius; 
for convenience we look at the case where 7=1. 

The proof in [5^ closely follows that of [TU] for conductivities a S W^'P{VI)^ 
p > 2. The matrix potential is, however, defined slightly differently, and 
since the potential in [2^ is not Hermitian, the approach in [20] is to consider 
the complex case as a perturbation from the real case provided the imaginary 
part of 7 is small. Define Q-yiz) and a matrix operator D by 

Thus we define 

Qi2{z) = -^9^1og 7(z) and Q2i{z) = -^9^ log 7(2), (6) 
and equivalently we can write 

= - Q2i(z) = — (7) 

Qi2(z) - - Vi(z)57(z) and Q^^^z) ^ -\r\z)d^[z) . (8) 



or 



Defining a vector 



in terms of the solution u to ([T|), one sees that 

\w j ' \w 
The uniqueness result in [2^ is 

Theorem 2.1. (Theorem 1.1 '201) Let Vl be an open hounded domain in 
with Lipschitz boundary. Let Uj and Cj, for j = 1,2 satisfy assumptions and 
\\(t\\w'^-°°{Q.)j Ikll < There exists a constant LOi) ~ a;o(/3, cq, 17) such 
that if 7j = (Tj + itucj for j = 1, 2 and uj < loq and if 

then 

(Ti = (72 and ei = £2- 
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2.1 CGO solutions 

Francini shows in [20j that for lo sufficiently small and 7 satisfying ([3]) and (|H) 
there exists a unique 2x2 matrix M{z, k) for fc G C satisfying 

Af (•, k)- I e LP(m2), for some p>2, (10) 

that is a solution to 

{Dk-Q^{z))M{z,k)^0, (11) 
where 1?^ is the matrix operator defined by 

DkM = DM - ik 

and "off" denotes the matrix consisting of only the off-diagonal entries of M. 
The system pT|) is equivalent to the following set of equations, included for the 
reader's convenience 



5.M11-Q12M21 = id,-lk)Mi2-Ql2M22 = 

{d, + ik)M2i - Q21M11 = d,M22-Q2iMi2 = 0. 

Thus, there exists a unique matrix 5'(z, k) defined by 



(12) 



vi/(z,fcj-Anz,fcj Q e--^j - l^e-'=M2i(z,fc) e--'=M22(z, fc) j ' ^^^^ 
that is a solution to 

(7? - Q^) * = 0, (14) 

or equivalently 

9.*11-Ql2*21 = a,*2i-Q2i*ii -0 (15) 

9^*12 - (3l2*22 = <9z*22 -Q21*12 = 0. 

These CGO solutions '^{z,k) are key functions in the reconstructions, but 
the proof in [50] does not provide a link from these functions to the DN data. 
A useful link can be established through exponentially growing solutions to the 
admittivity equation ([l}. For 7 — 1 with compact support, equation ([l} can 
be studied on all of R^, and introducing the complex parameter k, two distinct 
exponentially growing solutions, which differ in their asymptotics, exist. We 
will denote these solutions by ui and U2 where ui ~ and U2 ~ ^~ik~ ^ 
sense that is made precise in Theorems 12.21 and 12. 3[ where the existence of such 
solutions is established. The proof will make use of the following lemma proved 
in the real case by Nachman [43]; the complex version shown here also holds 
and was used in [3D]. The lemma is also true if dz is replaced by dz- 

Lemma 2.1. Let 1 < s < 2 and - ^ - - h- 

r s 2 

1. If the complex function v G L^(M^), then there exists a unique complex 
function u G L'^iJS?) such that [dz + ik)u = v. 
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2. If the complex function v G L''(R2) and d^v 6 L^(R2), k eC\ {0}, then 
there exists a unique complex function u S W^'^'{^) such that {dz+ik)u = 

V. 

3. If the complex function v G L''(M^) and dzV e L%R^), fc e C \ {0}, then 
there exists a unique complex function u G W^'^{^) such that {dz~ik)u = 

V. 



The following lemma will also be used in the proofs of Theorems 12.21 and \ 

Lemma 2.2. Foruj sufficiently small andj satisfying ^ and the following 
identities hold: 

d,{-,{zy'/Hln{z,k) - 1) = {d,+ik){^{z)-'/^M2i{z,k)) (16) 
dzhiz)-'/^M22{z,k) - 1) = {dz-ik)ijiz)-'/^M,2iz,k)). (17) 
Proof. By the product rule, 

4 {j{zy^^^Mn{z,k) - l) - 4 (tW"'/') Mn{z,k)+j{z)-^/^dMiiz,k)) 

= ^{z)-'/^Q2iiz)Mniz, k) + i{zr^'^Qi2{z)M2i{z, k) 

= 7(^)"'^'(9z +ifc)M2i(z, fc) + 7(^)"'/'Qi2(2)M2i(z, fc). 

The second and third equalities utilized ^ and (IT^ . respectively. 
We also have 

(a,+ifc) (7(z)-i/2^/2i(z,fc)) =9z (7(z)-i/'M2i(z,fc)) +zfc7(z)-i/2^/2i(z,fc) 

= a, (7(^)"'/') M2i(z, fc) + 7(^)"'/' 9z(M2i(z, fc)) 

+ ifc7(z)-i/2M2i(z,fc) 
= 7(2)"'/'Qi2(^)M2i(^, fc) + -i{z)-^/\d, +lk)M2l{z, fc). 

This establishes ([TC)) . 

Similarly, using and ([T^ . 

9, (7(z)-1/2M22(^, fc) - l) = (7(2)"'/') M22(Z, fc) + 7(2)"'/' a.(A/22(z, fc)) 

- 7(^)"'/'Ql2(^)M22(^, fc) + 7(^)"'/'Q2l(^)A/i2(z, fc) 

= 7(^)"'/'(4 -^fc)Ml2(z, fc) + 7(^)^'/'Q2i(^)Mi2(^, fc). 

We also have 

(4 -ik) (7(z)-i/2a/i2(z, fc)) = a, (^j{z)-^/Hh2{z, fc)) - ifc7(^)"'/'Mi2(z, fc) 

= 4 (7(^)"'''') Mi2(z, fc) + liz)-'^^ d,{M,2{z, fc)) 

-ifc7(z)-i/2Mi2(z,fc) 
= 7 W/'Q2i(^)A/i2(z, fc) + 7(^)"'/'(4 -~tk)Mi2{z, fc). 

This establishes (H?]). □ 
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Knudsen establishes the existence of exponentially growing solutions to the 
conductivity equation in the context of the inverse conductivity problem in |32| . 
The proofs of their existence for the admittivity equation and the associated 
boundary integral equations are in the same spirit as [32] . 

Theorem 2.2. Letj{z) G W^'P{fl), withp > 2 such that a ande satisfy ^ and 
and let j{z) — 1 have compact support in W^'P{fl). Then for all k E C \ {0} 
there exists a unique solution 



ui{z,k) = 6^'=^ 



— +wi{z,k) 
ik 



(18) 



to the admittivity equation in such that wi{-,k) £ W^''''{M.'^), 2 < r < oo. 
Moreover, the following equalities hold: 



id. +ik) 



-ikz ^ 



1 



e "'""ui(z, k) 

ik 



e *''^Mi(a;, k) 

ik 



^-i-^'^{z)Mii{z,k) ~ 1 
^-t-^'\z)M2i(z,k), 



and 



e ^'^''ui{x,k) 



<C 1 



1 



(19) 
(20) 

(21) 



for some constant C. 



Theorem 2.3. Let 7(2) satisfy the hypotheses of Theorem \2.2l Then for all 
fc e C \ {0} there exists a unique solution 



U2{z, k) = e 



-ikz 



— +W2{z,k) 
ik 



(22) 



to the admittivity equation in with W2{-,k) G W^'^{M.'^), 2 < r < 00. More- 
over, the following equalities hold: 



(9z -ik) 



e^^''u2{z,k)-^ 

IK 



ikz 



U2{z,k) 



ik 



= J-'/\z)M22iz,k)-l 

= 7-i/2(z)Afi2(z,fc), 



and 



e U2{z, k) 



1 



< C 1 



1 



(23) 
(24) 

(25) 



for some constant C. 

We will prove Theorem 12.21 the proof of Theorem 12.31 is analogous. 
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Proof. Assume u is a solution of the admittivity equation of the form (ITS)) , and 
let {v,w) = 7^/^(9z u, i9z u) be the corresponding solution to {D — Q-y)'^ = 0. 
Define the complex function v via v {z) = 7(z)-i/2Mii(z,fc) - 1. We will first 
show there exists a unique complex function wi S W^''^{M.^), where r > 2 such 
that {dz +ik)w = u, for fc G C \ {0}. Let us rewrite v as follows: 



v{z) = 7(z)-i/2[Afn(^,fc)-l 



7(z) 



-1/2 



1 



Let r > 2 and 1 < s < 2 with 



: i — i. We know by Theorem 4.1 of 
that there exists a constant C > depending on /3, ctq and p such that 



sup||Mii(z,fc) - 1| 



< C for every r > 2, and that 7(2) — 1 has com- 



pact support in T4^^'''(R^). It follows that v e L''(M^), and by Minkowski's 
Inequality 



< a 



\Hz)h,- = li^r'^' [^hi{z, fc) - 1] + [liz)-'^' - 1 

where Cr,^ depends on r and the bounds on a and e. 
From ®, 

= (4 7(^)^'/')Mn(^, fc) + 7(^)-'/'(4 Afii(^, fc)) 

= j{z)-'/^Q2iiz)Mniz,k) + ^{z)-'/^Qi2{z)M2i{z,k) 

= j{z)-^^^Q2i{z)[Mniz,k)^l] 

= l{zr^'^Ql2{z)M2l{z, fc) + l{z)-^'^Q2l{z) 

We know that -f{z)-^/'^Q2i{z) G ^"(M^) with I < a <p since (5i2(z) has 
compact support. It follows that 7(2)^^/^(321(2) £ ^''(M^) n L2(M2-)_ gy ^j^g 
generalized Holder's inequality and the fact that ||Mii(z,fc) — is bounded 
with 7 = 7 + 5, we have d.2.v{z) £ L*(R^) and ||4 i'||^3(][j2) ^ -f^r,7, where Kr^-y 
depends only on r and the bounds on cr and e. Thus, by Lemma \TJ\ (2), there 
exists a unique solution wi{z,k) G iy^'''(M^) such that 



(9z +ifc)wi(2, fc) = 7(z)-i/2Mii(z, fc) - 1. 

We have by P^ . 

a(7(z)-i/2^-fii(z,fc) - 1) = (a, +zfc) (7(z)-i/2Af2i(z,fc) 
Taking 4 of both sides of ([26| and using (|27|) , 

4 (a, +ifc) u;i(z, fc) = 4 (7(z)-i/2Mn(z, fc) - l) 
^{dz+tk) [j{z)-'/Hl2i{z,k) 
Using the fact d{d + ik) = (d + ik)d, it follows that 

id, +ik) (4 wi{z, fc) - 7(2)-i/2M2i(z, fc)) - 0. 



(26) 



(27) 



(28) 



(29) 
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Since d^wi{z,k) - -/{z)-'^/'^M2iiz,k) e L''{R'^), by Lemma [O (1), we must 
have 

dz wi{z, k)=j{zr'/^M2iiz,k). (30) 



We now define 



ui{z, k) — e 



ikz 



Wl{z, k) + 



(31) 



then by ^ 

{dz +ik) (^e-'^'ui{x, k)-^^^ (5z +ik) wi{z, k) = j-^^^{z)Mn{z, k) - 1 



which proves pJ)) . and by ([30 



[e-'^^'u^ix, k)-^^ 4 w^{z, k) = 7-1/2 (^)Af2i(z, 



which proves ((20)) . 

The norm estimate given by (j2ip fohows by Minkowski's Inequality, the 
constant C depends on r, the bound on 7 — 1, and the bounds on a and e. □ 

Remark: Note that from (H 



(32) 



7-i/2Mn(z, fc) - 1 ^ (a. +ifc) (^e-^'=^wi(z, fc) - ^ 



5z(e-'''"ui) + iA:e-"^Mi(z, k) - 1 



-2^2: 



dzUi{z,k) - 1, 



and from (l20 



7-i/2M2i(z,fc) (^e-^'=^ui(z,fc) - 1^ 

= u{z, k) 4 (e-*'^") + e-*'^" 4 ui(^, fc) 



(33) 



-ikz 



dzUi{z,k). 



Thus, we can equivalently rewrite (ITOl) and (1^ . respectively, as 

7i/2a,ui(z,fc) = e''=^Mn(z,fc) = ^'n(z,fc) 
7i/2(z)a,ui(z,fc) = e^'=^Af2i(^,fc) = *2i(^,fc). 



In a similar manner, we can rewrite ([23]) and ([24]), respectively, as 

7i/2(z) a, ^2(2, fc) = e-^'=^"M22(z, fc) ~ 2e-''=^>i/2 
-1^/\z)d,U2{z,k) = e-^'=^"Afi2(z,fc) = vI;,2(z,fc). 



(34) 
(35) 



(36) 
(37) 



Useful boundary integral equations for the traces of ui and U2 can be derived 
under the additional assumption that 7 e W^^^ and ui,W2 G W'^'P , p > 1. The 
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following proposition shows a relationship between the exponentially growing 
solutions ips{z,k) (when they exist) to the Schrodinger equation 

{-^ + qs{z))i;s{z,k)=Q (38) 

and the CGO solutions ui and U2 to ([1]). The solution tps to (1551) . where qs is 
complex, is asymptotic to e^^^ in the sense that 

ws = e-^'=^^s(-, k)~leLPr\ 

where h = h — \ and 1 < n < 2. The question of the existence of a unique 
solution to (|38l) is addressed for real 7 in [43], where it is shown to exist if 
and (roughly) only if qs = ^^1/2 ■ The solutions "05 will be used to derive the 
boundary integral equations for ui and U2, but not in the direct reconstruction 
algorithm. 

Lemma 2.3. Let 7(2) = a{z) + iLxje{z) G VF^'P(f2), with p > 2 such that a and 
e satisfy ^ and and let 7(0) — 1 have compact support in W^'P(fl). Let 
Ml be the exponentially growing solution to the admittivity equation as given in 
Theorem \2.2[ and let ips &e the exponentially growing solution to the Schrodinger 
equation psp . when it exists. Then 



ifcui(z,fc) = 7-i/2(z)V's(2,fc). (39) 

Proof. Note that 

ikui{z, k) = e^'^^il + ikwi{z, k)) 



e 



li^y^'^ [7(2)'/' + l{zY'^ikwi(z, k) 
^^7(z)^i/2 (1 + [7(2)1/2 _ ^^j ^ 7(z)i/2ifcwi(z, fc)) 



(40) 



satisfies the admittivity equation with 7(2) -'^/^ — l + ^i{zY/'^ikwi{z^ k) G W'^'^{Vt) 
for r > 2. We also know that when it exists 

7(2)-i/V5(^, k) = e^'=^7(^)"'/'(l + ws{z, k)), 

is also a solution to the admittivity equation with ws{z, k) e ^^^'^(R^). Hence, 
these exponentially growing solutions must be equal. □ 

We will recall some terminology arising from [43j before establishing bound- 
ary integral equations involving the exponentially growing solutions. The single 
layer potential Sk is a boundary integral operator defined by 

Skf{z)= [ f{OGk{z~OdsiO, (41) 
Jdn 

where Gk is the Faddeev Green's function defined by 

Gfc(z) := e*^gfc(z), -AG,. = <5, (42) 
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and 



(2^)2 



-d^, (-A - 4ikS)gk = S, 



(43) 



e(e + 2fc) 

for A; G C \ {0}. In the real- valued case 7 = cr, the trace of the function ijjs{-, k) 
on dfl satisfies the integral equation [33] 



^s{-, k)\dn = e"=^|an - 5fc(A, - Ai)^s(-, k), 



(44) 



where G C \ {0} , A^- is the voltage-to-current map when ft contains the 
conductivity distribution a, and Ai is the Dirichlct-to-Neumann map of the 
homogeneous conductivity 1. The equation is a Fredholm equation of the 
second kind and uniquely solvable in H^^^{dQ) for any k £ C \ {0} . 
The boundary integral equations for ui and U2 are similar to (I44p . 

Theorem 2.4. Let 7 G Vl^^'^(ri) /or p > 1 and suppose 7 = 1 m a neighborhood 
of dfl. Suppose a and e satisfy ^ and and let 7(2) — 1 have compact 
support in Then for any nonexceptional fc G C \ {0}, the trace of the 

exponentially growing solution ui{-,k) on dQ is the unique solution to 

ui(z,fc) = -Le'*^^ -S'fc(A^ - Ai)ui(-,fc), zedn. (45) 

Proof Let i = i - i, where 1 < r < 2 and p > 2. Let {7„}„gN C W^''"{n) be 
a sequence converging to 7 G W^''P{ft). Then by the Sobolev Embedding Theo- 
rem, {7„}„gpj C W^''^{fl). Let 4'n be the exponentially growing solutions to the 

Schrodinger equation with potential 7n^^^A7y^, and u" be the exponentially 
growing solutions to the admittivity equation with admittivity 7„. Then for 
each 71 G N, the complex 7 version of ((44|) holds for nonexceptional k G C \ {0} 



V'n(-,^)|aa 



e''laf2-^fc(A^„-Ai)V'„(-,A:), 



where 7„ = 1 in the neighborhood of dH.. 

It follows by ((39|) that for each complex number fc 7^ 0, and for each n G N 



In 



-1/2 



ik 



^„ = k) ui{-, fc) in H^^^idil). 



(46) 



We claim that for each ' 



satisfies (gS]). To see this, by 



ik 



5fc(A^-Ai)«"(.,fc) 



an 



1 

ik 
1 



Akz 



ikz 



5fc(A^-Ai)- 



,,-1/2 



an 



ifc 



-'0n(-,fc) 



5fc(A^ 



Ai)— V'„(-,fc) 
ik 



1pn{z, fc) 



9n 



7i 



-1/2 



ik 

u'\z,k), 



1pn{z, fc) 



an 



(47) 



11 



where we used the fact that 7^ = 1 in a neighborhood of dil. Thus, m" satisfies 
(|45)) for each n £ N. 

We know by Theorem 3.1 of [20] that M{z,k) depends continuously on 7. 
From pS)) . we can conclude that 

5fe(A^„ - Ai)w"(-, k) ^ Sk{A^ - Ai)ui(-, k). (48) 

Thus, by (gB]), (gT]), and (gH]), we have that mi(-, fc)lan satisfies (gS]). The unique- 
ness oi ui{-, k)\gfi follows by Theorem 12.21 □ 

An analogous theorem holds for 1*2: 

Theorem 2.5. Let 7 e W'^'P{fl) for p > 1 and suppose 7 = 1 m a neighborhood 
of on. Suppose a and e satisfy ([3]) o-nd (g]), and let 7(2:) — 1 have compact 
support in VF^'^(r2). Then for any nonexceptional k £ C \ {0}, the trace of the 
exponentially growing solution M2(-,A;) on dfl is the unique solution to 

U2{z,k) = ^e-''''' ~SkiAy~Ai)u2{;k), z e dn. (49) 
2.2 The Scattering Transform Matrix 

The scattering transform S-y{k) of the matrix potential is defined in |20] by 

TT Jk^ \-Q2iiz)e{z,k)Mii{z,k) J 1^ 

(50) 

where e(z, k) = exp{i(2:fc + zfc)}. Thus we are only concerned with computing 
the off-diagonal entries of 5'^, which we will denote by Si2{k) and S2i{k). 

Boundary integral formulas for the off-diagonal entries of S^{k) in ([50]) can 
be computed by integration by parts as follows 

S,2{k) - - [ Qi2{z)e-''^''^22iz,k)dfi{z) 

"*^"4*i2(2,fc) d^l{z) 



Jn 

^ I e~''''-^i2{z, k){v^ + iV2) dS{z) --Id, e-*'=^*i2(z, k) dfi{z) 
271- Jon 

i 

2^ 



e-'^'"*i2(z, k){vi + iy2) dS{z) (51) 

an 



and similarly, 



27r 



S2i{k) e»''*2i(2, fc)(;^i - 1^2) dS{z), (52) 



where v ^ 1^1 + ii'2 denotes the outward normal to the boundary 9f2. 
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Theorem 2.6. The trace of the exponentially growing solutions ^'12(2:,^) and 
\E'2i(2;, k) for k G C \ {0} can be determined by 



*12(2,fc) 

"^21(2, k) 



ao 47r(z - C) 



47r(z - C) 



[A^-Ai]u2(C,fc)d5(C) (53) 
[A^-Ai]ui(C,fc)d5(C), (54) 



where ui and U2 are calculated via equations (|45p and (I49p respectively. 

Proof. We use the relations in ((35)) and p7l) to obtain boundary integral equa- 
tions for 5*21 and '^12 for z G dfl from Equations (|45|) and (|49p . respectively. 
Let us begin with ^'12: 



^12(2, fc) = -/'^^{z)d,U2iz,k) 



1 



-ikz 



ik 



Gfc(z-C) [A^-Ai]«2(C,fc) dS{C) 



Y/2(z) / 9,[Gfe(z-C)][A^-Ai]«2(C,fc)rf^(C). (55) 



Similarly, 

«'2l(2,fc) 



7i/2(z)a,wi(z,fc) 



- ^,1/2 



ik 



an 



Gk{z-0 [A^-Ai]7.i(C,fc) d5(C) 



= -l^'^z) 4[Gfe(z-C)][A^"Ai]^/i(C,fc)d5(C). 



an 



From HSl 



a 1 e{z,-k) 
■^9k(z) = — tkgk(z) 

ox Attz 'ilTZ 



^y9k{z) 



1 



Aniz 



e{z, —k) 
Aniz 



kgk{z), 



we arrive at 



and similiarly 



dzgk{z - C) = - 



e(z - C, -fc) 
47r(0-C) ' 



dz 9k(z - 



1 



- ikgk{z - C). 



(56) 

(57) 
(58) 

(59) 
(60) 



47r(z - C) 

Since Gk{z) — e*'^^.gfc(z), we use ([59)) and (jSU)) to compute the and 9^ deriva- 
tives of Gk{z — C), obtaining 

9zGfe(z-C) - 9z fe^'(^-^)gfc(z-C) 



gife(z-C) 

'47r(z--C)^ 



(61) 
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and similarly, 









_47r(z - C). 



Thus, 



(62) 
(63) 



Substituting the representations for the and 9z-derivatives of G/j(z — C), given 
in (pT|) and (|^ . back into the equations for \l/i2 and \E'2i, given in (I55|) and 
([SS]) respectively, proves the theorem. □ 

2.3 From 5(A;) to M 

The dependence of M on the complex parameter k is related to the scattering 
transform through the following dk system. 

Theorem 2.7 (Theorem 4.1 [2IT). Let a and e satisfy ^ and (g]) and let 

M he the unique solution to (llip satisfying (jlOp . T/ie map fc — M{-,k) is 
differ entiahle as a map into L^_p, and satisfies the equation 



dkM{z,k) = M{z,k) 



e{z,k) 
e{z,~k) 



S^{k), 



(64) 



wher 



S^{k) = - 



e{z, —k) 






-e(z, k) 



{Q^M)''ff{z,k)dniz). (65) 



Moreover, for every p > 2, 



sup||M(z,-) -/||^p(jj2) < K2, 

k 

where K2 depends on (3, CTq; and p. 

Notice that Equation ([M)) can be written as the following two systems of 
equations: 

^kM^^{z,k) = A'h2{z,k)e{z,-k)S2i{k) 
9feMi2(z,fc)= Mii(z,fc)e(z,fc)5i2(fc) ' 



and 



dkM2i{z,k) = M22{z,k)e{z,-k)S2i{k) 
dkM22{z,k) = M2i{z~k)e{z~k)Si2{k) 



(66) 
(67) 



included for the reader's convenience. 
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2.4 From M to 7 

Theorem 2.8 (Theorem 6.2 For any p> 0, 



Q^iz)^ lira ^i{Bp{0))-^ DkM{z,k)dn{k). (68) 

'=0^°° J {k:\k-k„\<p} 

This provides a reconstruction formula for the entries of Q-y, and one can 
recover 7 from Q12 = — ^ log(7) or Q21 — —^d^ log(7). However, this formula 
is computationally impractical as it requires a large k limit of integrals involving 
Bz and dz derivatives of M{z, k). 

We have derived computationally advantageous formulas for recovering the 
entries of that only require knowledge of the CGO solutions at /c = 0. 
Theorem 12.91 provides this direct relation between the CGO solutions M(z,0) 
(from the dk equation (|64p ) and the matrix potential Q^{z)^ eliminating the 
large k limit required in equation (j68p above. 

Theorem 2.9. The entries of the potential matrix Q-y{z) defined in ([5]) can be 

calculated using only knowledge of the CGO solutions M{z, 0) via 

4m+(q^,z,o) 
^"^'^ ^ M_(g,,z,o) 

^"^'^ - M+(Q.,.,o) • 

M+(0^,z,fc) = Mii{Q^,z,k) + e{z,-k)Mi2{Qy,z,k) (71) 
M^{Q^,z,k) = M22{Q^,z,k) + e{z,k)M2i{Qj,z,k). (72) 

Proof. We follow an idea similar to that in [5 and define 

M+{Q^,z,k) = Afii(g^,z,fc) + e(z,-fc)A/i2(g^,2;,fc) (73) 
M_{Q^,z,k) = M22{Q^,z,k) + e{z,k)M2i{Qj,z,k). (74) 

Note that M+ and M_ are only dependent on the (^7 matrix, not —Q^ as is 
required in [5]. Therefore, 

dzM+{Q.y,z,k) = Qi2(z)e(z,-fc) [Af_(g7,z,fc) 

+ i(fc-fc)Mi2(Q7,z,fc)] 
dzM_{Q^,z,k) = Q2i(2)e(z,fc)M+(Q^,z,fc), 



so that 



4m+(Q7,z,o) = gi2(z)M_(g7,z,o) 

5zA/_(Q7,z,0) = g2i(z)M+(g7,z,0). 



□ 
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One can then reconstruct the log of the admittivity 7 from either Q12 or 
Q21 by inverting the or operator respectively, and exponentiate to recover 
7 explicitly 



log(7(2)) = 



Q2i(C) 



d/.(C). (75) 



2.5 The steps of the algorithm 

We now have all the necessary steps for a direct reconstruction algorithm: 

1. Compute the exponentially growing solutions ui{z, k) and U2{z, k) to the 
admittivity equation from the boundary integral formulas (1451) and p9|) 



ui{z,k)\on = — — 
ik 



dn Jdn 



Gfe(z-C)(A^-AiW(C,fc)d5(C) 



U2{z, fc)|sn 



-ikz 



ik 



an Jdn 



2. Compute the off diagonal entries of the CGO solution ^'(z, k) for z G dfl 
from the boundary integral formulas ([55]) and ([Ml) 



*12(2,fc) 
*2l(2,fc) 



an 47r(z - C) 

gifc(2-C) 



an 



47r(z-C) 



[A^-Ai]u2(C,fc) dSiO 
[A^-Ai]wi(C,fc) d5(C). 



3. Compute the off-diagonal entries of the scattering matrix Sj{k) from ([51 
and ([HI) 



^12 W = 7TZ / e-*'=^^'l2(z,fc)(^^l+^^^2)d5(z) 
/an 



^21 (fc) 



i 

'2^ 



e'^'^'2iiz,k){iyi-ii^2)dS{z). 



on 



4. Solve the dk equation ([64|) for the matrix M(z, k) 

dkM{z,k) = M{z,k) 

5. Reconstruct from Theorem 12.91 and use ([75]) to compute 7. 



e{z,k) \ o 
e(z,-fc) P^W- 
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3 Numerical Implementation 



In this section, we describe the implementation of the algorithm. Greater detail 
of the numerical methods and validations of the computations for admittivity 
distributions with twice differentiable real and imaginary parts can be found 
in |25| . where the solution to the forward problem (|lip is computed and used 
to validate formulas (15^ an ([M]). as well as computations of the scattering 
transform. In this work, we consider examples with discontinuities at the organ 
boundaries. 

3.1 Computation of the DN map 

An approximation to the DN map was computed by simulating voltage data by 
the finite element method (FEM), and then computing a matrix approximation 
to the map by computing the inner product of the applied currents with the 
voltages. This approximation to the DN map has been discussed, for example, 
in [27l [ITl |39] . It can be formed analogously in the complex case. 

Gaussian white noise was added independently to the real and imaginary 
parts of the simulated voltages for each current pattern by adding a random 
vector of amplitude 77 > multiplied by the maximum voltage value for that 
current pattern and real or imaginary component to the computed voltages. We 
consider noise levels 77 = and ry = 0.0001, which corresponds to 0.01% noise, 
the published level of the ACT 3 system [18] , which applies the trigonometric 
current patterns used in the simulations here. 

3.2 Computation of the CGO solutions and S^{k) 

The CGO solutions on the boundary of were computed for each fc in a grid 
[—K,K]'^ in the complex plane. The choice of K, which serves as a cut-off 
frequency, was determined by the behavior of the scattering transforms 5i2 and 
5*21. As in |39| for the D-bar algorithm for conductivity reconstructions, the 
cutoff frequency K has a regularizing effect, and was chosen here empirically to 
balance smoothing and numerical error. We do not address the selection of K 
by more sophisticated means in this work. 

3.2.1 Computation of ui and U2 

A boundary integral equation of the form (|45|) was solved in [17^ and [39* . In 
this work, as in |17j . we employ an approximation to the Faddeev's Green's 
function Gk that allows for very rapid computation of ui and U2 from (|45p and 
(|49p respectively. Namely, Gk is approximated by the fundamental solution for 
the Laplacian 

Go{z) = — log|2:|. 
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Denoting the solutions to (HSI) . (|49p by ul and Uj, respectively, the convolution 
integrals 



Ui{z,k)\gn = 



ik 

^ — ikz 



Go(z-C)(A^-AiK(C,fc)d5(C) 

on Jdn 



Go(z-C)(A^-Ai)uO(C,fc)d5(C) 
dn Jan 



were computed for z = Z£, the center of the ^th electrode, via Simpson's rule, 
and Go was set to when ( = zi. 

3.2.2 Computation of 'I'i2 and ^21 

The boundary integral formulas ([53|) and (f54| for \E'i2 and ^'21, respectively, 
require knowledge of [A-^ ~ Ai]uj{C,k) for j ~ 1,2, with ( 6 9r2, and k € 
C \ {0}. These values are already computed during the evaluation of ui and 
U2 via (|45)) and ([49]). Therefore, we merely recall those values and approximate 
the boundary integral using a finite sum. One should note that Gq{z — ^), 
Gk{z — (), and Gk{z — are all undefined for z = (. We removed these 
points in the computation by setting their values to zero. 

3.2.3 Computation of the scattering transform 

The off-diagonal entries of the scattering transform matrix, namely Si2{k) and 
S2i{k), were computed inside the square [—KjK]"^ (with fc = not included 
since the formulas for the CGO solutions do not hold for fc = 0). We compute 
Si2{k) and 5*21 (fc) using a finite sum approximation to (jSip and (j52p : 

L 



Si2{k) « |^e-*^-^^^'i2(z,,fc)z, 



=1 

L 



S2iik) « -f ^e*'"'vl/2i(z^,fc)z7, 



L 

1=1 

where z/ denotes the coordinate of the equally spaced electrode around d^l 
(in this case the unit circle). 

3.3 Solution of the system of D-bar equations 

The two systems of dk equations ([55)1 and ([57)) can be written as the convolutions 

1= Mn(z,fc)-^*(Mi2(z,fc)e(z,-fc)52i(fc)) , . 

0= Mi2(z,fc)-i*(Mn(z,fc)e(z,fc)5i2(fc)) ' 

and 

1= M22{z,k)-^*{M2i{z,k)e{z,k)Si2{k)) f^^. 
0= A/2i(z,fc)-^*(M22(z,fc)e(z,-fc)52i(fc)) • 
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A numerical solver for equations of the form 



afet;(fc) = T{k)v{k). 

was developed in [35 for the inverse conductivity problem. The solver is based 
on the fast method by Vainikko [36] that uses FFT's for solving integral equa- 
tions with weakly singular kernels. 

In this work, we must solve the systems of equations ([75]) and ([77)1 rather 
than a single equation. Furthermore, the unknowns M(z, k) are not conjugated, 
but instead the argument k is conjugated. To address this, we interpolated 
the scattering data S*^, computed above in Section |3.2.3[ to a new fc-grid that 
includes the origin fc = at the center and has an odd number of grid points in 
both the horizonal and vertical directions. We solve the systems ((7S)) and (|77p 
on this new /c-grid using appropriate flip operations to ensure that we access 
the correct entries in the matrix corresponding to Af(z, k). 

To perform the convolution we used Fourier transforms as follows: 



1 

— * 

nk 



and similarly 
1 

— * 
tt/c 



{Ah2{z,~k)e{z,-k)S2i{k)) 

hi IFFT (fFT ■ FFT {M^^iz, k)e{z, -k)S2i{k)) ) , 



{Mn{z,k)e(z,~k)Si2{k)) 



= hi IFFT (fFT (^-^^ • FFT (Mn(z, k)e{z, k)Si2{k))^ , 

where ft.^ is the step size of the uniform /c-grid of size 129 x 129, and • denotes 
componentwise multiplication. We used GMRES to solve the resulting linear 
systems for each value of z in a grid of 128 equally spaced points between [- 
1.2,1.2] in both the x and y directions and computed M{z,k) for all \z\ < 1.2. 
The step size in z was /i^ = 24/1270 « 0.019. 



3.4 Computation of the admittivity 

The admittivity is computed by solving first for Q21 from (|70|) (note that equiv- 
alently one could use Q12 from and then solving ([75t for log(7) in the 

Fourier domain using FFT's. The functions M+ and M_ in equations (|7T|| and 
([7^ were evaluated using the entries of M{z,0) recovered when solving the dk 
equation (see Section 1531 above). We used centered finite differences (with a 
step size of h^ = 24/1270 w 0.019) to evaluate the and dz derivatives of A/+ 
and M_ respectively. We then performed componentwise division to compute 
Q12 and Q21 for \z\ < 1.2. Finally, the computed log(7) was exponentiated to 
recover 7 inside the unit disk. 



19 



4 Numerical Results 



We consider several test problems simulating a simplified cross-section of a hu- 
man torso. In each example, the admittivity is given by 7 = a + ie. That 
is, the imaginary component includes the temporal angular frequency uj. Since 
this is a known value, there is no loss of generality in representing 7 this way 
in the simulations. The complete electrode model (CEM), originally described 
in |14j . was implemented in the FEM in order to solve the forward problem. 
The CEM takes into account both the shunting effect of the electrodes and the 
contact impedances between the electrodes and tissue. In our computations, 
was chosen to be a disk of radius O.ISto, and the FEM computations were per- 
formed on a mesh with 4538 triangular elements and 32 equispaced electrodes 
0.029to X 0.024m placed on the boundary. The effective contact impedance was 
chosen to be 2: = 0.0057f2m^ on all electrodes in our simulations. The current 
amplitude was chosen to be C = 2mA, and the applied current patterns are the 
trigonometric patterns 



Ccos{j9i), 1<£<L, l<j<f 

Csin((f 1<£<L, f + l<j<L-l, 



where 9£ = ^jf-, je^j is the area of the electrode, le is the current on the 
electrode, and L denotes the total number of electrodes. As in [371 HZ], the 
currents were normalized to have £^-norm of 1, and the voltages were normalized 
accordingly. Also, the DN map was scaled to represent data collected on the 
unit disk using the relation A-y^i = rAj^r, where the second subscript represents 
the radius of the disk. 

Where indicated, we added 0.01% Gaussian relative noise to the simulated 
voltages as follows. Denote the (complex-valued) vector of computed voltage 
for the j-th current pattern by , let r/ = 0.0001 denote the noise level, and N 
a Gaussian random vector (generated by the remdn commmand in MATLAB) 
that is unique for each use of the notation N. Denoting the noisy data by 
we then have — Ke{V^) + i lni{V^) where 

Re{V^) ^ Re{V^)+r]ma.x\Re{V^)\N 
Im{V^) = Im{V^)+7]max\lm{V^)\N. 

We solve the boundary integral equations (|45|) and p9)) for the traces of the 
CGO solutions ui and U2 for k G [—K, K]'^, with K varying for each test prob- 
lem in this work. The solution M{z,k), to the dk equation ([M]) . is computed 
in parallel by the method described in Section 13.31 The low-pass filtering by 
taking k € [~K, K]'^ results in smooth functions Mjp, j,p = 1, 2, which are dif- 
ferentiated by centered finite differences to recover Q21, as described in Section 
13.41 The admittivity 7 was then computed by (TTS]) . 

Define the dynamic range of the conductivity, and likewise the permittivity, 

by 

max CT^^' — min cr*^^) ^ , , 
: 100%, (79) 
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where the maximum and minimum values are taken on the computational grid 
for the reconstruction and cr^^^ denotes the reconstructed conductivity a that 
was computed using a scattering transform computed on the truncated k grid. 

4.1 Example 1 

The first test problem is an ideahzed cross-section of a chest with a background 
admittivity of l+Oi. We do not include units or frequency in these examples, 
since our purpose is to demonstrate that the equations in this paper lead to a fea- 
sible reconstruction algorithm for complex admittivities. Reconstructions from 
more realistic admittivity distributions or experimental data are the topic of 
future work. Figure 1 shows the values of the admittivity in the simulated heart 
and lungs. Noise- free reconstructions with the scattering transform computed 
on a 128 X 128 grid for k G [—5, 5]^ are found in Figure 2. The reconstruction 
has a maximum conductivity and permittivity value of 1.1547 -f 0.1422i, occur- 
ring in the heart region and a minimum of 0.8497 — 0.0345i, occurring in the 
lung region, resulting in a dynamic range of 76.2428% for the conductivity and 
47.4027% for the permittivity when the negative permittivity value is set to 0. 
Although this decreases the dynamic range, we set the permittivity to when it 
takes on a negative value in any pixel, since physically the permittivity cannot 
be less than 0. The reconstruction has the attributes of good spatial resolution 
and good uniformity in the reconstruction of the background and its value. 



0.8 + 0.1 




1 +0 i > 



Figure 1: The test problem in Example 1. 



4.2 Example 2 

This second example was chosen with conductivity values the same as in Exam- 
ple 1, but with permittivity values in which the "lungs" match the permittivity 
of the background. This is motivated by the fact that at some frequencies, 
physiological features may match that of the surrounding tissue in the conduc- 
tivity or permittivity component. This example, purely for illustration, mimics 
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Figure 2: Reconstruction from noise-free data for Example 1 with the real part 
of 7 (conductivity) on the left, and the imaginary part (permittivity) on the 
right. The cut-off frequency was K — 5. The dynamic range is 76% for the 
conductivity, and 47% for the permittivity. 



that phenomenon. The admittivity values can be found in Figure 3. Noise-free 
reconstructions with the scattering transform computed on a 128 x 128 grid for 
k £ [—5, 5]^ are found in Figure 4. The maximum value of the conductivity 
and permittivity occur in the heart region, 1.1517 -I- 0.1583z, and the minimum 
value of the conductivity and permittivity is 0.8477 — 0.0222i. In this exam- 
ple, the dynamic range is 75.9850% for the conductivity and 52.7516% for the 
permittivity when the negative permittivity value is set to 0. Again the spatial 
resolution is quite good, and the background is quite homogeneous, although 
some small artifacts are present in both the real and imaginary parts. 

4.3 Example 3 

Example 3 is an admittivity distribution of slightly higher contrast, and a non- 
unitary background admittivity of 70 — 0.8 + 0.3*. See Figure 5 for a plot of 
the phantom with admittivity values for the regions. Due to the non-unitary 
background, the problem was scaled, as was done, for example, in ^17 ] I27 j . by 
defining a scaled admittivity 7 = 7/70 to have a unitary value in the neighbor- 
hood of the boundary and scaling the DN map by defining A.,^ = 7oA^, solving 
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o.a + i 




+ i 



Figure 3: The test problem in Example 2. Notice that in this case, the permit- 
tivity of the lungs matches the permittivity of the background, and so only the 
heart should be visible in the imaginary component of the reconstruction. 




0.85 0,9 1 1.1 1.15 Q.05 0.1 0.15 



Figure 4: Reconstruction from noise-free data for Example 2 with the real part 
of 7 (conductivity) on the left, and the imaginary part (permittivity) on the 
right. The cut-off frequency was K = 5. The dynamic range is 76% for the 
conductivity, and 53% for the permittivity. 

the scaled problem, and rescaling the reconstructed admittivity. The scattering 
data for the noise-free reconstruction was computed on a 128 x 128 grid for 
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k S [—4.7,4.7]. Noisy data was computed as described above in the beginning 
of this section, and the scattering data was computed on a 128 x 128 grid for 
k S [—4.3,4.3]. The reconstructions are found in Figure 6. The maximum and 
minimum values are given in Table 1. In this example, for the noise- free recon- 
struction, the dynamic range is 60.7503% for the conductivity and 60.9469% for 
the permittivity. Again the spatial resolution is quite good, and there is little 
degradation in the image and the reconstructed values in the presence of noise. 
The scattering transform did begin to blow up for noisy data, requiring us to 
choose a smaller square for the computation of the scattering transform than in 
the noise-free case. Choosing a square with if = 4.3 results in a dynamic range 
of 60.3647% for the conductivity and 55.9110% for the permittivity. Choosing a 
square with K = 4.5 results in a dynamic range of 66.8521% for the conductivity 
and 57.8504% for the permittivity with a small loss of resolution compared to 
the reconstruction with a K = 4.3. For the choice K = 4.5, this loss was very 
small, but for K = 4.7 it already resulted in a significant amount of distortion 
of organ boundaries. A thorough study of the effects of the choice of K and its 
method of selection is not included in this paper. 





Admitivity of 
test problem 


Reconstruction from 
noise-free data 


Reconstruction from 
noisy data 


heart 
lungs 


1.2 + 0.6 i 
0.5 -1- O.li 


1.0072+0.4573 i (max) 
0.5819 + 0.1526 i (min) 


1.0061 + 0.4389 i (max) 
0.5835 + 0.1593 i (min) 



Table 1: Maximum and minumum values in Example 3 with the non- unitary 
background were found in the appropriate organ region. The table indicates 
these values of the admittivity in the appropriate region. 



0,6 + 0.1 I 




+ 0. 1 i 



Q8 + 0.3i > 



Figure 5: The test problem in Example 3. In this case, the background admit- 
tivity is 0.8 -|- 0.3i, rather than 1 -I- Oi as in Examples 1 and 2. 



24 



5 Conclusions 



A new direct method is presented for the reconstruction of a complex conduc- 
tivity. This method has the attribiites of being fully nonlinear, parallclizablc, 
and the direct reconstruction does not reqiiire a high frequency limit. It was 
demonstrated on numerically simulated data representing a cross-section of a 
human chest with discontimious organ boundaries that the method yields re- 
constructions with good spatial resolution and dynamic range on non-noisy and 
noisy data. 
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Figure 6: Top row: Reconstruction from noise-free data for Example 3. The 
cut-off frequency was K — 4.7. The dynamic range is 61% for the conductivity, 
and 61% for the permittivity. Bottom row: Reconstruction from data with 
0.01% added noise. The cut-off frequency was K = 4.3. The dynamic range is 
60% for the conductivity, and 56% for the permittivity. 
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